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^ ! ABSTRACT 

The Munich Near-IR Cluster Survey (MUNICS) is a wide-area, medium-deep, photometric 
survey selected in the K' band. It covers an area of roughly one square degree in the K' 
and J near-IR pass-bands. The survey area consists of 16 6' x 6' fields targeted at QSOs with 
redshifts 0.5 < z < 2 and 7 28' x 13' stripes targeted at 'random' high Galactic latitude fields, 
ff^ , Ten of the QSO fields were additionally imaged in R and /, and 0.6 deg^ of the randomly 

^SJ ' selected fields were also imaged in the V, R, and / bands. The resulting object catalogues 

I were strictly selected in K', having a limiting magnitude (50 per cent completeness) of K' ^ 

19.5 mag and J ^ 21 mag, sufficiently deep to detect passively evolving systems up to 
a redshift of z < 1.5 and luminosity of 0.5L*. The optical data reach a depth of roughly 
, i? 23.5 mag. The project's main scientific aims are the identification of galaxy clusters 

pi i' at redshifts around unity and the selection of a large sample of field early-type galaxies at 

< z < 1.5 for evolutionary studies. In this paper - the first in a series - we describe the 
survey's concept, the selection of the survey fields, the near-IR and optical imaging and data 
reduction, object extraction, and the construction of photometric catalogues. Finally, we show 
^ I the J — K' vs. K' colour-magnitude diagramme and the R — J vs. J — K', V — I vs. J~K', 

and V — I vs. V — R colour-colour diagrammes for MUNICS objects, together with stellar 
population-synthesis models for different star-formation histories, and conclude that the data 
] set presented is suitable for extracting a catalogue of massive field galaxies in the redshift 

j_j ■ range 0.5 ^ z ^ 1.5 for evolutionary studies and follow-up observations. 
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mology: observations 



1 INTRODUCTION 

Directly observing the evolution of individual galaxies with time 
is, unfortunately, not possible. Therefore we must rely on investi- 
gating the statistical properties of the whole galaxy population as 
a function of redshift, trying to draw conclusions from ensemble 
properties on the evolution of typical members of these ensembles, 
and thus facing difficulties, like, for example, discriminating be- 
tween luminosity evolution and number density evolution. 

Much work has been invested in this field, resulting in a lot of 
progress in the last decade which has seen many imaging and red- 
shift surveys being undertaken using different selection techniques 
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in wave-bands from the UV to the sub-mm. These surveys have a 
wide range of scientific applications, from the detection of high- 
redshift galaxy clusters to the study of the evolution of 'normal' 
field galaxies. 

The earlier opticall y and near-IR selected redshift and im aging 
surveys, among others ^roadhure^^llis^^^han^^^8^)j_Col- 
less et al. (1990), |LiUy. Cowie, & Gardner (1991D, Glazebrook 
et al. (1994), and |Cowie et al. (1994| ), lai d the path to the lan dmark 
CFRS (Canada-France Redshift Survey; Lilly et al. 1995a), an I- 



band selected redshift survey mapping the evolution of the galaxy 
population out to z ~ 1. Also many 'pencil-beam' surveys have 
been carried out, t he most prominent being the Hubble Deep Field 



North and South ( [Williams et al. 1996^ [Williams et al. 200C| ) and 



their ground-based imaging and spectroscopic follow-ups, allow- 
ing us a first glimpse at the galaxy population at 2.5 ^ z ^ 4.5. 
The inability to determine redshifts spectroscopically for all ob- 
jects where multi-band imaging data is available (because of lim- 
ited telescope resources, either in observing time or in collecting 
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area), caused ph otometric redshift deter mination techniques to gain 
attention again^^^^^63^^^^85|^Fernandez-Soto, Lanzetta, 
& Yahil 1999r^enfte^00^!'Thisrand the wide-field imagers be- 
coming available in the optical and also in the near-infrared wave- 
length regime, made multi-band imaging surveys a very promising 
option for further studies in galaxy evolution. 

Selection in in a single pass-band introduces different (and 
sometimes subtle) selection effects, a well-known fact which need 
not necessarily be considered at the disadvantage of the resulting 
object database, as long as the selection function is well-understood 
and under control. These selection effects can be used deliberately 
for probing different galaxy populations and different aspects of 
their evolution. While selection in blue pass-bands is used to study 
star forming sources, selection in the near-IR is predominantly sen- 
sitive to the light of old stellar populations. Near-IR fc-corrections 
are small even at redshifts above unity and insensitiv e to the spec- 
tral type of the observed objects ( jCowie et al. 1994 ) and to short- 
lived bursts of star fo rmation, as has been pointed out by K auff- 
mann & Chariot (1998). Thus near-IR selected surveys are thought 
to be much less biased with respect to the mix of spectral types 
compared to optically selected surveys. Furthermore, the uncertain- 
ties resulting from inhomogeneous dust absorption are minimal in 
the near-IR. It has therefore been concluded that near-IR selection 
is a feasible attempt at a selection in stellar mass (R ix & Rieke 
1993; [Brinchmann & Ellis 2000| ). 

The Munich Near-Infrared Cluster Survey (MUNICS) is an 
attempt at closing the gap between previously undertaken infrared- 
selec ted deep pencil-beam surveys (Gardner, C owie, & Wainscoat 
1993 ; iMcLeod et al. 1995|; |Cowie et al. 1994 D jorgovski et al. 
1995; |Williams et al. 1996|;^aracco et al. 1997|) and relative ly shal- 
low wide-area surve^s^^MobasheTj^^hamleSj^^^l^^ 
brook et al. 1994; pardner et al. 1997[ ), simultaneously profiting 
from the advantages of near-infrared selection. 

MUNICS is a wide-area, medium-deep, photometric survey 
selected in the K' band. One part of the surveyed fields was centred 
on known quasars, while the rest was randomly selected at high 
Galactic latitudes. It covers an area of roughly one square degree 
in the K' and J bands with optical follow-up imaging in the I, R, 
and V bands for a large fraction of the total surveyed area. 

The resulting object catalogues are strictly selected in K' with 
a limiting magnitude of K' ~ 19.5 mag and J ~ 21 mag, suffi- 
ciently deep to detect passively evolving systems up to a redshift 
of 2 ^ 1.5 and luminosity of 0.5L* (see Fig. [l]). The optical data 
reach a depth of roughly R ~ 23.5 mag. 

This paper is laid out as follows. Section]^ deals with the sur- 
vey's scientific aims and concept. The required sensitivity in terms 
of limiting magnitudes, the selection of the survey fields as well 
as the photometric system adopted are described. In Section ^ we 
give an overview of the observations and discuss the reduction of 
the near-IR and optical imaging data. Section ^ contains the data 
analysis. The methods for object detection, photometry, and object 
classification are discussed, and a first analysis of the survey's com- 
pleteness as well as number counts in five colours are presented. 



2 SURVEY CONCEPT AND LAYOUT 
2.1 Scientific aims 

The project's main scientific aims are the following. First, to iden- 
tify clusters of galaxies at high redshift by detecting their lumi- 
nous early-type galaxy population. As has been shown in the last 
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Figure 1. Stellar population synthesis models in the J — K' vs. K' plane 
for different star-formation histories, a SFR(i) oc S(t) burst at z = 4, 
and 3 exponential star formation rates SFR(t) oc exp(t/r) with t = 1, 
3, and 10 Gyr, setting in at 2 = 4. The models are normalised to have 
a luminosity of Lg at the present epoch, according to the type-dependent 
luminosity function of the Virgo cluster (see text) as given by Sandage, 
Binggeli, & Tammann (1985). The thin dotted line indicates the limiting 
depth (50 per cent completeness) of the MUNICS data in J and K'. The 
adopted cosmology is Hq = 65, Hq = 0.3, Q/^ = 0.7. The SSP models 
used in the synthesis are from Maraston (1998). Additionally, an SSP model 
by Bruzual & Chariot (1995) is shown for comparison as a thin solid line. 



years, the early-type gala xy population in clusters is well in place 



at redshifts of at least 0.8 ([Stanford, Eisenhardt, & Dickinson 1995 
Stanford et al. 1997^ [Stanford, ^jsenh^^^^DkMnso^^9^rde 



Propris et al. 1999). Given the small fc-corrections in the K band, 
this makes selection in the near-IR a promising approach to detect 
clusters at redshifts around unity, complementing selection in other 
optical bands. Clusters of galaxies allow to find large numbers of 
massive galaxies at higher redshift and thus represent unique lab- 
oratories to study the evolution of galaxies in high-density regions 
as a function of redshift, and in contrast to the evolution of similar 
galaxies in the field. Furthermore, the evolution of the number den- 
sity of clusters is a promising test of cosmological models, depend- 
ing sensitively on the density parameter Qp (Eke, Cole , & Frenk 
1996 ; [Bahcall, Fan, & Cen 1997|; [Bahcall & Fan 1998t Eke et al. 
1998). While the number of clusters known at redshifts z > 0.5 
is steadily increasing (mostly due to X-ray selection), samples se- 
lected uniformly in the optical and near-IR wavelength ranges are 
still deficient. 

Cluster detection at high redshifts is strongly biased towards 
the most massive systems, mainly because of lack of detection sen- 
sitivity for lower mass systems. Finding also less massive systems 
is important when reasoning about hierarchical galaxy formation 
models, since the galaxies in the densest environments formed ear- 
lier, so by looking only at the most dense environments one is ef- 
fectively pushing the epoch of collapse, merging, and star forma- 
tion out to higher redshifts and further away from the observational 
window. Therefore we decided to centre a subset of the MUNICS 
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fields on known quasars hoping to increase the chance of detecting 
clusters in their environment. 

Secondly, a statistically well-defined sample of the early-type 
galaxy population in the field can be constructed from our cata- 
logues, which will be used to study the evolutionary history of such 
objects in the redshift range < 2 < 1 by means of the iC-band se- 
lected luminosity function, the luminosity density at near-infrared 
wavelengths, and the two-point correlation function. Again, K- 
band selection offers unique opportunities due to the close connec- 
tion betw een near-IR luminosity and stellar mass (B rinchmann & 
Ellis 2000), and thus allows direct assessment of the predictions of 
hierarchical galaxy formation theories. 

Thirdly, the natureof^xtremel^^e^_object^^ 
Rieke, & Rieke 1988; |Hu & Ridgway 1994 ) will be examined. 



EROs, usually defined in terms of R — K greater than or ap- 
proximately equal to 5 at moderately faint if -band magnitudes of 
K > 18, are thought to be either high- redshift early-type g alaxies 
or heavily extincted starburst galaxies (|Cimatti et al. 1999|; Small 
et al. 1999), the relative contribution of the two sub-populations 
being yet highly uncertain. Due to the small areas of the surveys 
available so far, even the surface de nsity of these objects is not reli- 
ably known (Thompson et al. 1999). Since they mostly are 7?-band 



'dropouts', having the possibility to detect such objects in the MU- 
NICS data in the / and J bands, together with the large field cov- 
ered, will enable us to gain valuable information on their nature. 



2.2 Limiting sensitivity 

Fig. |l| shows stellar population synthesis models in the J — K' vs. 
K' plane for different star-formation histories, a SFR(f) oc 5{t) 
Simple Stellar Population (SSP), and 3 exponential star formation 
rates SFR(t) oc exp(t/r) with t — 1,3, and 10 Gyr. The onset of 
star formation occurs at z = 4 in all models. The populations are 
normalised to have a luminosity of L*g at the present epoch, accord- 
ing to the ty pe-dependent luminosity function of the Virgo cluster 
as given by [Sandage, Binggeli, & Tammann (1985 ). The adopted 
values for Mg are Mg = -21.5 for the SSP model (elliptical/SO 
galaxy), Mg = -20.5 for r = 1 (Sa-Sb spiral), Mg = -19.5 
for r = 3 (Sc), and = -17.5 for r = 10 (Sd and later). The 
cosmology is Ho = 65 km s^^ Mpc~ \ ^Ip = 0.3, Qa = 0.7. 
The SSP models are taken from 



Maraston (\99i). The dis- 



tinguishing feature of that synthesis method is the adaption of the 
fuel consumption theorem to evaluate the energetics of the post 
main-sequence evolutionary phases. The models used here have 
solar metallicity and age ranging from 30 Myr to 15 Gyr. The Ini- 
tial Mass Function (IMF) is a power law *(A/) oc M"'^+''' with 
Salpeter exponent x = 1.35 down to a lower mass limit of O.IA^q. 
The optical and infrared colours predicted by these SSP models are 
calibrated against Milky Way and Magellanic Cloud globular clus- 
ters and com£ared^o_simi^^ 

ton (1 998). Fig. |l|also shows an SSP model by B ruzual & Chariot 
(1993) using the 1995 version of their code, with solar metallicity 
and a Salpeter IMF. The models evolve similarly up to redshifts 
of '--^ 1. The differences in colour are likely due to the different 
treatment of the post main sequence stages and are discussed in 
Maraston (199^ . 

Following the predictions of these models, the limiting magni- 
tudes in the near-IR wave-bands have been chosen to be 19.5 mag 
in K' and 21.0 mag in J, such that early-type objects having lu- 
minosities of ^ 0.5L* at the present epoch can be detected in K' 
virtually at any redshift, and in J up to a redshift of 2: ^ 1.5, as- 
suming passive evolution. 
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Figure 2. Layout and nomenclature of one of the MUNICS 'mosaic fields'. 
The size of the stripe is 13' X 28' , covered by 2 X 4 pointings in the near-IR 
and two pointings in the optical with the circular field of view of CAFOS 
(see text). The IR image is divided for technical reasons into two 2x2 
mosaics. Each optical frame and IR mosaic frame are denoted by the name 
of the stripe, here S2, followed by the subfields they cover, giving S2 fl. . . f4 
and S2 f5. . . f8 for this stripe. 



This is in agreement with the findings of the CFRS, which has 
shown that, while the luminosity function of the population of blue 
field galaxies shows significant signs of evolution in the redshift 
range 0.2 < 2 < 1 - explainable by brightening or increase in 
space density - the redder part of the population (roughly redder 
than She) shows no signs of evolution of its luminosity function 
in the same redshift range ( Lilly et al. 1995b ). The latter is inter- 
preted in terms of brightening of the individual galaxies through 
passive evolution counterbalanced by negative density evolution, 
such that the luminosity function of the early-type population ef- 
fectively does not evolve. 



2.3 Field selection 

The MUNICS survey consists of two sets of near-IR target fields, 
one set of single camera pointings having an effective field of view 
of 6' X 6' pointed towards quasars, and a second set of 28' x 13' 
fields constructed from mosaics of pointings targeted at random 
high Galactic latitude fields. This second set of fields was selected 
to contain no bright stars, nearby bright galaxies, and known nearby 
clusters of galaxies, and furthermore, to have low Galactic redden- 
ing (which is all together difficult to accomplish together with the 
prerequisite of having no bright star within the field, given our field 
size). 

A total of 16 fields targeted towards quasars with redshifts 
0.5 < 2 < 2 were observed. These fields will be referred to as 
'quasar fields' hereafter, labelled Ql. . . Q16. The quasars were se - 
lected from the seventh edition of the Veron-Cetty & Veron (1996) 



catalogue. The selection criteria were B < 19.0 mag, 0° < Dec 
< 65°, 0*" < RA < 18*", and 0.5 < 2 < 2. Six of these quasars are 
not detected in the radio bands of the catalogue (6 cm and 1 1 cm) 
and are therefore considered radio quiet. The remaining 10 are ra- 
dio loud. 

A second set of 7 fields was targeted at high Galactic latitude 
'empty' fields, i.e. free of bright stars {V < 17 mag) and known 
nearby extragalactic objects. These fields will be called 'mosaic 
fields' hereafter, for they are mosaiced images in the near-IR. They 
are labelled S 1 . . . S7. Each such field is laid out as a stripe of 4 x 2 
IR pointings, yielding an area of 28' x 13'. For technical reasons, 
namely that four near-IR pointings can be completed in K' and J 



4 A^. Drory et al. 




0.8 1.0 



Figure 3. Relative transmission of MUNICS V, R, I, J, and K' filter (solid 
curves) curves including quantum efficiency of the CCD (V,R, and /) 
and Rockwell HAWAII near-IR array (J and K'), as well as the atmo- 
spheric transmission in the near-IR. Relative transmission curves of stan- 
dard Johnson-Kron-Cousins V, R, and / filters, as well as J and K' are 
shown for comparison (dotted curves). 
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during a single night as well as image size and efficiency of optical 
follow-up observations (see below), each such stripe is divided into 
two 2x2 mosaics of single IR frames, denoted fl~f4 and f5-f8. 
This particular geometry was chosen for efficiency since such a 2 x 
2 mosaic of IR frames suits the circular field of view of the optical 
imager we used (CAFOS, having a diameter of roughly 16 arcmin; 
see Section 3.2). For clarity, Fig. ^ shows a sketch of the geometry 
and nomenclature of the mosaic fields. 

Finally, Table |l| lists the coordinates, available pass-bands, 
seeing, and Galactic foreground extinction of the observed mosaic 
and quasar fields. 



2.4 Photometric system 

The MUNICS imaging observations were carried out partly using 
non-standard filters, or imperfect realisations of standard filters (see 
Section^ below). Since the colours of objects in the MUNICS cata- 
logues extend to much redder colours than any available photomet- 
ric standard stars, we decided to work in the MUNICS instrumen- 
tal photometric system and not to transform magnitudes into the 
standard Johnson-Kron-Cousins system. Linear transformation to 
the Johnson-Kron-Cousins system would have caused magnitude 
errors up to 1 mag, because the true transformations are highly 
non-linear, especially for red objects. The MUNICS photometric 
zero-points are in the Vega system. 

Note that, since comparison of the object's colours with spec- 
tral synthesis models is intended (e.g. for deriving photometric red- 
shifts or discussing the nature of EROs), it is important that the ob- 
served colours and the synthetic colours are consistent with respect 
to the filter set. 

Accurate measurements of the transmission curves of the glass 
filters and quantum efficiencies of the detectors were obtained and 
applied in all subsequent synthetic photometry. The filter curves are 
shown in Fig. ^ 

Such accurate knowledge of the filter system allows a reli- 
able calibration of different bands via colour-colour diagrammes 
of stars which are compared with synthetic stellar sequences ob- 
tained from the convolution of SEDs from stellar libraries with the 



3f^ 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — 



o „o o c 




+ ■ ++ ■ X * A 



0.0 



0.5 



1.0 



R-I 



1.5 



2.0 



Figure 4. Comparison oi R — J vs. V — R (upper panel) and J — K vs. 
R — I (lower panel) colours derived by convolving stellar SEDs (see text) 
with the MUNICS filter curves (open circles) with a sample of bright stars 
(R. < 21) identified in the MUNICS fields S2 f l-f4 (crosses), S4 f l-f4 (x), 
S6 f5-f8 (squares), and S7 f5-f8 (fi-iangles). 



transmission curves. The absolute photometric zero-points can then 
be derived from a single photometric observation in one band only. 

Fig. ^ shows a comparison of _R— J vs. V — R and J — K\s. 
R — I colours derived by convolving stellar SEDs with the MU- 
NICS filter curves with a sample of stars detected in the MUNICS 
mosaic fields. The stellar SEDs used for computing the synthetic 
colours are taken from the Bruzua l-Persson-Gunn-Stryker spec- 
tral library ( |Gunn & Stryker 1983 ; Strecker, Erickson, & Wi tte- 
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bom 1979), covering spectral types 05 to M8. The agreement be- 
tween the synthetic photometry and the data along the stellar loci in 
colour-colour space demonstrates that the constructed filter curves 
match the actual ones and that the photometric zero-points are mu- 
tually consistent in the optical and the near-IR regime (see also 
Sect. 3.2 and Sect. 3.1). It is worth noting that we have no cool gi- 



ants or supergiants in the MUNICS sample (as those would occupy 
the redder sequence in J—K at i?— J ^ 1.5) and no stars of earlier 
type than roughly late F to G. 



3 OBSERVATIONS AND DATA REDUCTION 

3.1 Infrared observations and data reduction 

The A'' -band an d J-band imaging was ob tained using the Omega- 
Prime camera (Bizenberger et al. 1998) at the prime focus of 



the Calar Alto 3.5-m telescope. Omega-Prime is equipped with a 
HAWAII 1024^ HgCdTe array. The image scale is 0.396 arcsec 
per pixel, resulting in a 6.75' x 6.75 ' field of view. The K' fil - 
ter (Ao = 2.12^i, AA = 0.35^i; see [Wainscoat & Cowie I992| ) 
was used because it significantly reduces the thermal background 
seen by the detector relative to the standard K filter, thus gaining 
sensitivity. Table ^ lists all observing runs undertaken to present 
date. 

The A'' -band data were observed using a dithering pattern 
consisting of 16 positions within an area of 30" x 30" laid out on a 
4x4 grid with 10" spacing between adjacent grid points. The data 
were recorded using a randomized sequence of these 16 positions. 

On each position 28 seconds of net exposure time were col- 
lected, divided into several shorter exposures as necessary depend- 
ing on the ambient temperature and thus the level of the thermal 
background. The length of the single exposures was always chosen 
such that non-linearity of the detector was negligible. 

This 16 position cycle was repeated 3 times yielding 48 frames 
and a total exposure time of 1344 s. The J-band images were 
observed using the same dithering pattern with longer integration 
times of 80 s on each position and therefore needed only one cycle, 
giving a total of 1280 s. 

The near-IR mosaic fields consist of four such Omega-Prime 
pointings arranged in a 2 x 2 configuration with 6' offset in each 
direction measured from field centre to field centre. Each mosaic 
then covers a total area of 162 square arc minutes, counting only 
the central area with the longest total exposure time and removing 
overlaps and borders due to the dithering pattern (see Fig. 

On photometric nights, standard stars from the UK IRT Faint 
IR Standard Stars catalogue (Casali & Hawarden 1992) were ob- 



served several times during the night at different air masses to 
determine the photometric zero point and the atmospheric extinc- 
tion coefficient. To increase the number of standard star measure- 
ments available for e ach night, the calib rations of further stars in 
the UKIRT fields by iHunt et al. (19981) were included. Night-to- 



night variations in the zero-point were typically less than O.I mag. 
Targets observed during non-photometric nights were re-observed 
(with shorter exposure time) at least once during photometric con- 
ditions to assure accurate photometric calibration. The typical for- 
mal uncertainties in the zero-points were 0.05 mag in A' and 
0.06 mag in J. The extinction coefficients were found to be sta- 
ble for all runs with typical values around 0.08 ± 0.025 mag per 
airmass in A' and 0.12 ± 0.02 mag per airmass in J. By compar- 
ison with synthetic photometry as explained in Sect. 2.4 we con- 
clude that additional systematic errors in the near-IR calibration 



as well as systematic offsets between the near-IR and the optical 
wave-bands cannot be larger than '--^ 0.1 mag. 

The data were reduced using standard image processing algo- 
rithms within IRAlQ. For each frame a sky frame was constructed 
from typically 6 to 12 (temporally) adjacent frames where bright 
objects and detector defects have been masked out, and which were 
scaled to have the same median counts. These frames were then 
median-combined using clipping to suppress fainter sources and 
otherwise deviant pixels to produce a sky frame. The sky frame 
was scaled to the median counts of each image before subtraction 
to account for variations of sky brightness on short time-scales. The 
sky-subtracted images were flat-fielded using dome flats to remove 
pixel-to-pixel fluctuations in quantum efficiency. The frames were 
then registered to high accuracy using the brightest ~ 10 objects 
and finally co-added, again using clipping to suppress highly de- 
viant pixels due to cosmic ray events and defective pixels on the 
array, after being scaled to airmass zero and to a common photo- 
metric zero-point. 

The 2x2 mosaic images were produced by registering the 
images using objects in the overlap regions, simultaneously cross 
checking the photometric calibration. Before combining, the im- 
ages were adjusted to have the same background counts computed 
from the mode of the pixel values in 'empty' sky regions of the 
images to correct for residual differences in sky brightness. Abso- 
lute astrometric calibration of the images is discussed in Sect. 3.3 
below. 



3.2 Optical observations and data reduction 

Optical imaging of the mosaic fields was performed at the Calar 
Alto 2.2-m telescope in the V, R, and / bands using the Calar Alto 
Faint Object Spectrograph (CAFOS) focal reducer in direct imag- 
ing mode. CAFOS was equipped with a SITe 2048^ CCD detector, 
yielding a resolution of 0.53 arcsec per pixel and a circular field 
of view (due to vignetting by optics) of 16' in diameter. The V- 
band filter used was a standard Johnson filter, the _R-band filter was 
an R2 filter (Ao = 0.648/x, AA — O.I68/1), slightly narrower and 
bluer than Kron-Cousins R. The /-band filter was an RG780 fil- 
ter with the red cutoff set by the CCD (see Fig. ^.) Total exposure 
times were 2700 s in V and /, and 1800 s in R, divided into several 
shorter exposures taken with offsets of ~ 15", depending on the 
presence of bright stars and on seeing conditions to avoid too many 
saturated objects. 

The quasar fields were imaged using the Imaging Grism In- 
strument (IGI) at the 2.7-m telescope of McDonald Observatory, 
using a 1024^ TK4 CCD (7' field of view) and Mould R and I in- 
terference filters. Exposure times 1800 s in i? and 2700 s in /, again 
divided into several shorter exposures. 

The optical CCD data were reduced in a fairly standard man- 
ner using IRAF, except for cosmic ray cleaning. The frames were 
bias/overscan corrected and then flat-fielded using a combination 
of dome flats and sky flats. The /-band frames showed consider- 
able fringing. Fringe images were created from the affected series 
of science exposures and occasionally also from twilight flats by 
medianing de-registered images after masking bright sources by 
hand as necessary. In some cases it was necessary to subtract a low 



IRAF is disti'ibuted by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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Field coordinates are given wit h respect to the image centres . 
QSO designations according to Veron-Cetty & Veron (199(^ . 



(1) Mosaic incomplete in the near-IR. 

(2) Near-IR data quality poor. 

(3) Good data quality in all five wave-bands. 

(4) Radio quiet QSO. 



order fit to the overall background in the science frames prior to 
construction of the fringe image to account for changes in the illu- 
mination pattern present in the images in the case where a bright 
star was close to the image border. The fringe images were then 
appropriately scaled and subtracted from the affected frames. 

Cosmic ray events were identified by searching for narrow lo- 
cal maxima in the image and fitting a bivariate rotated Gaussian 
to each maximum. A locally deviant pixel is then replaced by the 
mean value of the surrounding pixels if the Gaussian obeys ap- 
propriate flux ratio and sharpness criteria (Gossl, Riffeser, & Fliri 
2001). Such a procedure is much more expensive in terms of com- 
puting time (roughly 10 CPU minutes per frame) compared to stan- 
dard median filtering techniques, but is much more reliable in find- 
ing cosmic ray events in the wings of objects and in cleaning long 
cosmic ray trails. 

The re-imaging system of CAFOS causes substantial radial 
distortion of the image which had to be dealt with before co-adding 
the offset images. Therefore the frames were rectified using the 
known distortion equation, a polynomial of fourth order in the dis- 
tance from the optical axis (K. Meisenheimer, private communica- 
tion). 

If necessary, variations in the background intensity across the 
frames caused by scattered light were fitted and subtracted in each 



individual frame. The images were then corrected for atmospheric 
extinction and scaled to a common photometric zero-point before 
finally being added using the positions of ~ 15 bright objects for 
determination of the offsets between the individual frames. 

Duri ng photometric nights, photometric sta ndard stars were 
observed (Landolt 1992; Christian et al. 1985) and programme 
fields with insecure calibrations were re-observed with short expo- 
sures. The run at the Wendelstein 0.8-m telescope was devoted to 
such re-calibration to have independent calibrations for the fields. 
For each field, a photometric zero point and the atmospheric extinc- 
tion were determined. No colour terms were fitted to the calibration 



data, as explained in Sect. 2.4 The typical formal uncertainties in 
the photometric calibration were ~ 0.08 mag in /, ~ 0.04 mag in 
R, and ~ 0.05 mag in V. The extinction coefficients were usually 
consistent with a Rayleigh atmosphere, with a few nights showing 
higher extinction, albeit within the variations typical for Calar Alto. 
As with the near-IR observations, consistency was checked against 
synthetic colours of stars, again finding no systematic offsets. We 
conclude that systematic errors in the optical photometric calibra- 
tion are again smaller than 0.1 mag in all filters. 
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Table 2. MUNICS observing runs 



Date 


Icl. 


Instrument 


Remarks 


1996 24-27.10 


CA35 


n' 


Quasar fields 


1997 15-19.5 


CA35 


n' 


Quasar fields 


1998 8-14.4 


CA35 


n' 




1998 12-17.5 


CA35 


n' 




1998 28.5-1.6 


CA22 


CAFOS 




1998 16-18.11 


McD27 


IGI 


Quasar fields 


1998 16-20.12 


CA22 


CAFOS 




1998 23-30.12 


CA35 


n' 




1999 18.3 


Wdst 


MONICA 


Calibrations 


1999 27.5-3.6 


CA35 


SI' 




1999 9-18.6 


CA22 


CAFOS 




2000 26-31.5 


CA35 


MOSCA 


Spectroscopy 


2000 27-28.5 


HET 


LRS 


Spectroscopy 


2000 16.7 


CA35 


n' 


(1) 


2000 20-22.11 


ESO-VLT 


FORSl/2 


Spectroscopy 


2000 24-28.11 


CA35 


MOSCA 


Spectroscopy 


2000 5.12 


CA22 


CAFOS 


(1) 


2000 17-18.12 


CA22 


CAFOS 


(1) 


2000 19.12 


CA35 


n' 


(1) 


2001 17-21.1 


CA35 


MOSCA 


Spectroscopy 


2001 11-13.2 


CA35 


n' 


(1) 



(1) Re-imaging of fields with poor data quality. 

CA22 and CA35 ai'e the 2.2-m telescope and the 3.5-m telescope of Calar 
Alto Observatory, respectively. McD27 is the 2.7-m telescope and HET 
is the Hobby-Eberly Telescope, both of McDonald Observatory, Austin. 
Texas. Wdst is the 0.8-m telescope of Wendelstein Observatory operated 
by the Universitats-Stemwarte Miinchen. 



3.3 Astrometry 

Astrometxic solutions were computed for all JC'-band images 
to translate pixel coordinates into celestial coordinates. For this 
purpose, astrometr ic standards from the USNO-SAl.O catalogue 
( Monet et al. 199( ) were selected in each frame. The celestial coor- 



dinates of these stars were matched against the pixel position using 
the IRAF task CCXYMATCH, and the plate solution was computed 
using CCMAP. The typical scatter is less than 0.4 arcsec rms. 

The V, R, I, and J images of each field were registered 
against the K'-bmd image by matching the positions of ~ 200 
bright homogeneously distributed objects in the frames and deter- 
mining the coordinate transform from the A'' -band system to each 
image in the other four pass-bands using the tasks XYXYMATCH 
and GEOMAP within IRAF. The scatter in the determined solutions 
is less than 0.1 pixels rms in the transformation from K' to J, and 
less than 0.2 pixels rms from K' to the optical frames. Note that 
the frames themselves are not transformed. We only determine ac- 
curate transformations and apply these later to the apertures in the 
photometry process. 



3.4 Spectroscopy 

A spectroscopic follow-up programme is currently being conducted 
at the Calar-Alto 3.5-m telescope, the Hobby-Eberly Telescope, 
and the VLT, aiming ultimately at a magnitude-limited redshift sur- 
vey of the A"'-band selected catalogue. The results of these obser- 
vations will be discussed in a future paper. In the mean time, the 
spectroscopic redshifts are used to calibrate photometric redshifts. 



4 DATA ANALYSIS 

The construction of photometric catalogues from the reduced im- 
ages will be discussed in this chapter. The individual steps in this 
process will be described in some detail. Thes e inc lude the detec- 
tion of objects in the A''-band images (see Sect. 4. 1 ), a first analysis 



of the survey's completeness (Sect, i.2), the photometry of objects 
in all filters (see Sect. [4.3| ), and the separation of stars and galaxies 
in the catalogue (see Sect. [4.4| ). The chapter will be concluded by 
a comparison of number counts for galaxies with previous studies 
(Sect. as a further consistency check on our data set. 



4.1 Object detection 

Object detecti on was per formed using the YODA source extrac- 
tion software ( ^rory 2001 ). This package was specifically designed 
to be used in multi-band imaging surveys, where the background 
noise is often inhomogeneous across the images - in mosaiced 
frames or in dithered images where the exposure time is a function 
of position - and where the frames do not share a conmion coordi- 
nate system and pixel scale, due to the use of multiple telescopes 
and imagers. The second point was considered a serious problem 
since re-sampling the images to a common coordinate system in- 
troduces quite considerable noise for faint sources. 

Sources are detected by requiring a minimum number Npi^ 
of consecutive pixels to lie above a certain threshold t expressed in 
units of the local rms a of the background noise. To foster detec- 
tion of faint sources, the images are convolved with a Gaussian of 
FWHM equal to the seeing in the image. The choice of the number 
of consecutive pixels Npix and the threshold t is somewhat a trade- 
off between limiting magnitude at some completeness fraction, say 
50 per cent, and the num ber of tolerable spurious detections per unit 
image area (Saha 1995). Note that, since we aim at purely A'-band 



selected catalogues for most of our applications, the presence or 
absence of a source in other wave-bands cannot be used for confir- 
mation or rejection of sources, so the expected number of spurious 
detections per unit image area is of great interest to us. 

To find reasonable values for Npix and t we performed simu- 
lations on the A'-band image of one of our mosaic fields (S6 f5- 
f8). The dependence of the 50 per cent completeness limit on Npix 
and t was determined by adding point sources to the A'' -band im- 
age and recording the fraction of the objects recovered by the de- 
tection software as a function of Npix and t. The number of false 
detections was determined by looking for positive detections in an 
inverted (multiplied by —1) version of the image, after convincing 
ourselves that the background noise was sufficiently well approxi- 
mated by a Gaussian. Fig. ^ shows the results of these tests. 

At the depth of our data we detect roughly 1000 objects per 
mosaic field. Accepting 1 per cent contamination by false de- 
tections, i.e. roughly 10 false objects per mosaic field we fixed 
the detection threshold at t = 3cr and the minimum number 
of consecutive pixels at 1.4 times the seeing disk area, Npix = 
1.47^(FWHM/2)^ (10 pixels at 1", 16 pixels at 1.5"seeing for 
the near-IR frames) and performed object detection using these pa- 
rameters on all A'-band images. 



4.2 Completeness 

To estimate the completeness limits of the MUNICS A'-selected 
catalogues, Monte-Carlo simulations were carried out to determine 
the detection completeness as a function of magnitude in each A'- 
band image. 
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Figure 5. The behaviour of the 50 per cent completeness hmit for point-hke sources and the number of spurious source detections as a function of the detection 
threshold t in units of the local background rms a and the required number Npi^ of consecutive pixels above the threshold in units of the seeing disk area 
7r(FWHM/2)^. The left panel shows the change in limiting magnitude at 50 per cent completeness as a function of Npi^ for detection thresholds of 2.0cr 
(solid line), 3. Oct (dashed line), and 4.00- (dotted line). The right panel shows the number of spurious sources integrated over all magnitudes per image (one 
mosaic field), again as a function of Npi^ and the detection threshold t. Line styles as in the left panel. 



Table 3. Completeness limits for the MUNICS mosaic fields 
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S6 fl 
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S6f5 
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Band 


90% 


50% 
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90% 


50% 


90% 


50% 


90% 


50% 
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19.05 


19.45 
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19.18 


19.42 


19.80 


18.26 


18.56 


19.53 


19.92 


19.23 


19.83 


J 
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20.68 
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I 
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22.35 
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R 
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22.94 


23.32 
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23.58 
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23.37 


22.88 


23.28 


V 


23.13 


23.51 
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23.94 


23.36 
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23.36 


23.68 


23.50 


23.89 
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23.28 



In these simulations, 250 artificial objects with a Moffat-type 
PSF having the same FWHM as stars in the frames were added to 
the MUNICS images. For this purpose the IRAF package ARTDATA 
was used. A constant distribution of apparent magnitudes of the 
artificial objects was applied. Then the object detection algorithm 
was run on these frames, and the number of re-detected objects was 
recorded as a function of magnitude. This procedure was repeated 
250 times in order to decrease statistical errors. 

Fig. ^ shows the results of these simulations for 6 of the mo- 
saic fields where imaging data in all five colours are available. The 
50 per cent and 90 per cent completeness limits of the i^'-band 
images are listed in Table ^ From these simulations we conclude 
that the mosaic fields comprise a reasonably homogeneous data set, 
with field-to-field variations in 50 per cent completeness of order 
~ 0.4 mag, with the exception of S6 fl-f4 which is considerably 
shallower. The quasar fields, which are not shown here, have com- 
pleteness levels in the same range as the shown mosaic fields. 

Simulations with extended objects having de Vaucouleurs and 
exponential surface brightness profiles have also been performed. 
In general they yield 50 per cent completeness limits between 0.5 
and I magnitude brighter than those determined for stellar sources, 
depending mostly on the adopted profile scale-length. The full re- 
sults of these simulations will be extensively discussed in a future 
paper in the context of surface-brightness selection effects. 



4.3 Photometry 

Photometry was performed in elliptical apertures the shape of 
which were determined from the first and second moments of the 
light distribution in the /(''-band image, as described in D rory 
(2001), and additionally in fixed size circular apertures of 5 and 
7 arc seconds diameter. To ensure measurement at equal physical 
scales in every pass-band, the individual frames were convolved to 
the same seeing FWHM, namely that of the image with the worst 
seeing in each field. The signal-to-noise ratio as a function of mag- 
nitude for the 5"apertures is shown in Fig. ^ At the 50 per cent 
completeness limit in the K' band (19.59), the signal-to-noise ra- 
tio is ~ 10. For such an object having an R — K' colour of 6, the 
signal-to-noise ratio in the R band is roughly 3. 

Aperture fluxes and magnitudes were computed for each ob- 
ject present in the K'-band catalogue irrespective of a detection 
in any other band. For this purpose the centroid coordinates of the 
sources detected in the A''-band images were transformed to the 
other frames using the full astrometric transformations as deter- 
mined in Sect. [3.3| . The shape of the apertures were transformed 
using only the linear terms of the transformation. 

The photometric accuracy for the 5"aperture magnitudes is 
roughly O.I mag at K' = 19 mag. This error estimate includes the 
effects of photon noise and uncertainty in background determina- 
tion and subtraction, but does not include (systematic) errors due to 
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Figure 6. Completeness fraction of point-like sources as a function of mag- 
nitude in the K'-hand mosaic fields as determined by adding Moffat-type 
objects to the images and recording the fraction of such objects that was 
recovered by the detection process. 
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Figure 8. Magnitude errors as a function of magnitude in the mosaic 
field S6 f5-f8 in K' , J, I, R, and V as measured in circular apertures of 
5"diameter for i<''-band selected objects. Aperture fluxes are measured in 
every pass-band for each object present in the X'-band catalogue irrespec- 
tive of a detection in any other band. 
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Figure 7. Signal-to-noise ratio in K' (left panel) and R (right panel) as a 
function of magnitude in circular apertures of 5"diameter for objects taken 
from the i^'-selected catalogue of the field S6 f5-f8. The signal-to-noise 
ratio is defined here as the signal-to-noise ratio of the aperture photometry, 
i.e. total (sky-subtracted) flux within the aperture divided by the total noise 
within the aperture, with contributions to the latter coming from Poisson 
fluctuations in the object as well as the background, and the error in the 
determination of the background. 



the photometric cahbration. Fig. ^ shows plots of the magnitude er- 
ror vs. object magnitude for one mosaic field in all five pass-bands. 



4.4 Star-galaxy separation 

Star-galaxy separation relies on YODA's image classification stage 
which is based on a Bayesian analysis of the probability that an 
object's light distribution is due to an unresolved (point-like) source 
by comparison with light distributions constructed from the image's 
PSF. YODA's classification parameters are calculated for all objects 
in the catalogue in all avai lable pass-ban ds. 

As demonstrated in Drory (2001), classification is reliable 
across wave-bands and imaging instruments, and stellar sources al- 



most do not scatter out of the stellar locus in parameter space (ex- 
cept in the presence of crowding). Rather, images of faint galaxies 
as they become smaller at larger distances, move onto the locus of 
point-like sources. 

Using the multi-pass-band information available in MUNICS 
allows us to push the limit of reliable classification by using 
for each object the classification information in those pass-bands 
where the signal-to-noise ratio is highest. Therefore, in the mosaic 
fields where 5 colours are available, we classify as stellar every 
source that is classified as stellar by YODA in the three pass-bands 
with highest signal-to-noise. In the quasar fields, where less colour 
information is available, we rely on the two images with highest 
signal-to-noise. 

As can be seen in Fig. [l^ objects classified as stars occupy 
the clearly defined stellar sequence in the R—J vs. J — K' colour- 
colour plane, with only very few objects classified as stellar having 
a. J — K' colour redder than ~ 1. These are either misclassified 
faint and compact galaxies or very late-type stars or brown dwarfs, 
the latter is a possibility for those objects having also red R~ J 
colour. The objects lying on the stellar sequence at i? — J <; 2 and 
which are classified as galaxies were found to be faint and barely 
resolved objects failing the classification as a star only due to their 
appearance in one filter. In many cases an obvious reason - like a 
second close object - could be identified. We conclude that most of 
these objects are, in fact, misclassified stars. The total fraction of 
point-like sources in the catalogues is ~ 10 per cent. 

We have also checked the results of the image-based classifi- 
cation against spectral classification for those objects where spec- 
troscopy was already available, namely 45 galaxies and 53 stars 
having R < 20.5. All these objects were correctly classified. 
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Figure 9. Number counts for galaxie.s in K' , J, I, R, and V from MU- 
NICS data (filled circles) and previous studies (open circles), as described 
in the text. The counts shown are average number counts from all available 
MUNICS data and have not been corrected for incompletness. Error bars 
indicate Poisson errors. 



4.5 Galactic extinction 

We use the Galactic reddening_jna£sjrovide^^^ 
Finkbeiner, & Davis (1998) using a value of Rv ~ 3.1 to cal- 
culate A\ = R\E{B — V) to correct the measured magnitudes for 
Galactic foreground extinction. The values of E{B — V) for our 
fields are given in Table 0. 



4.6 Galaxy number counts 

In Table ^ we present number counts of galaxies in the MUNICS 
mosaic fields in all five filters K', J, I, R, and V. These counts 
are also shown in fig. ^, together with a compilation of number 
counts from the literature. Object catalogues were generated inde- 
pendently for each pass-band for this purpose, and star-galaxy sep- 
aration is based on the PSF classification as described above, using 
only single pass-band information. The data were not transformed 
into the standard magnitude system for this comaprison. Complete- 
ness corrections were not applied to these galaxy number counts, 
but corrected counts will be presented in the context of a more de- 
tailed completeness analysis in a future paper. The counts are av- 
erage counts from all the available MUNICS mosaic fields, with 
field-to-field variations in the number counts being on the level of 
0.1 dex. The errors given in Table ^ only include Poissonian errors. 

Th e number counts are compared to the following literature 
values. 3ardner, Cowie, & Wainscoat 1993; |Cowie et al. 199% 
Glazebrook et al. 1994|; [Djorgovski et al. 199^ ; [Gardner et al. 1996 
for t he K band, [Saracco et al. 1999|; Teplitz, McLean, & Malkan 



1999 for the J band, [Tyson 1988|; [Lilly, Cowie, & Gardner 1991 
Casertano et al. 1995|; [Gardner et al. 1996[ ; [Williams et al. 1996 
Huang, Cowie, & Luppino 1998[for the I band. Couch & Newell 



Koo 19861; [Stevenson, Shanks. & Pong 198^; |Yee & Green 1987 



1984; [Hall & Mackay 1984[; [Infante, Pritchet, & Quintana 1986 



Tyson 1988|; [lones et al. 199 1[; Metcalfe et al. 199 ![; [Picard 1991 



Couch, lurcevic, & Boyle 1993; [Steidel & Hamilton 1993[; Driver 
et al. 1994; [Metcalfe, Pong, & Shanks 1995| ; [Metcalfe et al. 1995 



1996 for the V band. 

We generally find good agreement with previously published 
number counts in all pass-bands, again as a consistency check con- 
firming the quality of our photometry. 



4.7 Colour distributions and objects at 2: > 1 

In Pig.|l^ we show the J~K' vs. K' colour-magnitude diagramme 
and the R-Jvs. J-K', V-I vs. J~K', and V - 1 vs. V-R 
colour-colour diagrammes for MUNICS data from three mosaic 
fields. The total number of objects shown is 2977, of which 286 are 
classified as point-like. These plots also contain the tracks defined 
by the stellar population synthesis models described in detail in 
Sect. [2.2[ . Briefly, the models are an SSP, and three exponential star 
formation histories with e-folding times of 1, 3, and 10 Gyr forming 
at 2; = 4. The models have been normalised such that they represent 
typical L* objects at z = 0, with L* chosen according to their 
'photometric' Hubble type. The cosmology adopted is again Ho = 
65,no = 0.3, J7a = 0.7. 

These models reasonably envelope the region in the colour- 
magnitude J~K' vs. K' plane occupied by the data, with the SSP 
model following the outline of the data points along the bright and 
red edge as might be expected since any further star formation or a 
later formation epoch would render the object bluer relative to the 
SSP 

It is also worth noting that the models constitute a continuous 
sequence with the duration of the star formation as the parameter 
in the R—Jws. J~ K' plane, closely following the SSP track up 
to a redshift of ~ 1, then rapidly turning bluer in i? — J while still 
getting redder in J — K'. A significant fraction of objects between 
the SSP and the 1 Gyr track is compatible with being well evolved 
objects at a redshift z ^ 1. How many objects exactly populate this 
region is an important question which will be addressed in a future 
paper. 

We finally conclude from these diagrammes that the quality 
of our data meets the requirements expressed in Sect. ^ and that we 
are in a position to construct a catalogue containing a large number 
of massive field galaxies in the redshift range 0.5 ^ z ^ 1.5 to 
study their evolution in detail. 



5 SUMMARY 

The Munich Near-IR Cluster Survey (MUNICS) is a wide-area, 
medium-deep, photometric survey selected in the K' band. It cov- 
ers an area of roughly one square degree in the K' and J near-IR 
pass-bands with additional complementary optical imaging in the 
V, R, and / bands. 

MUNICS has been undertaken to study the evolution of both 
field galaxies and galaxy clusters out to redshifts around unity, and 
to investigate the nature of extremely red objects and their connec- 
tion to the population of massive field spheroidal galaxies. 

The survey area consists of 16 6' x 6' fields targeted at QSOs 
with redshifts 0.5 < z < 2 and 7 28' x 13' stripes targeted at 
'random' high Galactic latitude fields. Ten of the QSO fields were 
additionally imaged in R and /, and 0.6 deg^ of the randomly se- 
lected fields were imaged in the V, R, and I bands. The result- 
ing object catalogues were strictly selected in K', having a limit- 
ing magnitude (50 per cent completeness) of K' ~ 19.5 mag and 
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Table 4. K', J, I, R, and V-band galaxy number counts for the MUNICS mosaic fields. The counts as a function magnitude and the error of the counts are 
given in logarithmic units. The values have not been corrected for incompleteness. The errors are Poissonian errors only. 
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J ~ 21 mag, sufficiently deep to detect passively evolving early- 
type systems up to a redshift of z ^ 1.5 and luminosity of 0.5L* . 
The optical data reach a depth of roughly R ~ 23.5 mag. The 
project's main scientific aims are the identification of galaxy clus- 
ters at redshifts around unity and the selection of a large sample of 
field early-t3T)e galaxies at < a < 1.5 for evolutionary studies. 

In this paper we describe the selection of survey fields as well 
as the observations and the reduction of the near-infrared and opti- 
cal data. We define our photometric system and show it to be inter- 
nally consistent by checking it against synthetic photometry of stars 
from stellar libraries. The construction of the if '-selected object 
catalogue is described in detail, particularly the choice of param- 
eters for object detection, which ensures completeness to as faint 
magnitudes as possible while keeping the rate of false detections 
in our strictly AT'-selected catalogue small. Photometry of the ob- 
jects in the catalogue is performed in elliptical apertures on frames 
convolved to the same PSF in all filters, in order to guarantee mea- 
surements of the flux in equal physical areas. Stars and galaxies in 
the fields are classified using a Bayesian analysis of the light distri- 
bution in the images. 

The quality of the survey data in terms of signal-to-noise ratio 
and limiting magnitude is discussed, with the completeness of the 
survey fields being characterised by Monte-Carlo simulations of 
point sources in the MUNICS frames. Also, galaxy number counts 
are presented in the five filters K', J, I, R, and V and compared to 
counts published by other authors. 

Finally, we show the J — K' vs. K' colour-magnitude dia- 
granmie and the J? — J vs. J — K' , y — 7 vs. J — K' , and V — I 
vs. V — R colour-colour diagrammes for MUNICS objects, to- 
gether with stellar population-synthesis models for different star- 
formation histories and conclude that the data set presented is suit- 
able for extracting a catalogue of massive field galaxies in the red- 



shift range 0.5 z ^ 1.5 for evolutionary studies and follow-up 
observations. 
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